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£T) [ On the lattice we study the gauge boson propagator of three dimensional compact 

QED in Landau gauge at zero and non-zero temperature. The non-perturbative 
qq ■ effects are taken into account by the generation of a mass, by an anomalous dimension 

and by the photon wave function renormalization. All these effects can be attributed 
to the monopoles: they are absent in the propagator of the singularity-free part of 
the gauge field. We assess carefully the Gribov copy problem for the propagator and 
the parameters emerging from the fits. 
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> : 

•'-j ■ I. INTRODUCTION 

Three-dimensional compact electrodynamics (cQED 3 ) has two essential features in com- 
mon with QCD, confinement 0] and chiral symmetry breaking 0. Although the physics 
behind it might be very different, monopole dynamics in three and in four dimensions, it is 
amusing to study certain nonperturbative aspects within a lower-dimensional model such 
as CQED3. Apart from its role as a toy model for QCD, the nonperturbative properties 
of CQED3 deserve attention in themselves because this model was shown to describe some 
features of Josephson junctions 0] and high-T c superconductors 0]. CQED3 was the first 
example, in space-time dimensions greater than two where it becomes a nontrivial problem, 
in which confinement of electrically charged particles was understood analytically It is 
the result of the dynamics of monopoles which emerge due to the compactness of the gauge 
field. Other common features of cQED 3 and QCD are the existence of a mass gap and of a 
confinement-deconfinement phase transition at some nonzero temperature. 

In two recent papers we have demonstrated how the deconfinement phase transition, 
occurring under the influence of compactification of one dimension, in 2 + 1 dimensions, 
finds its explanation from the monopole point of view [Hf and why the deconfinement phase 
transition is independent on the strength of the external fields 
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In a recent letter |7[ we answered the question what effect the confinement property has 
on the gauge boson propagator in this theory and what of the propagator is changing at 
the deconfinement temperature. The effects are twofold: first, an anomalous dimension 
appears which modifies the momentum dependence, and second, a mass is generated which 
can be well understood in terms of Polyakov's theory ^ . In Ref . Q we also have used the 
unique possibility in lattice simulations to remove the monopole degrees of freedom from 
the quantum gauge field in order to show that all nontrivial effects reside exclusively in the 
singular fields of the monopoles. 

Note that a non-trivial anomalous dimension of the gauge boson propagator may also 
appear due to dynamical matter fields [if. The sign (in our notation) of the anomalous 
dimension induced in this way is different and this leads to the binding of the monopoles 
into dipole pairs and, consequently, to the disappearance of the confinement at certain 
distances between the test particles. This picture was confirmed for the compact model in 
the presence of the dynamical matter fields 0. 

In the present paper we want to extend the analysis of Ref. [7[ in various respects. Since 
the propagator is studied in the Landau gauge, first of all we have investigated more carefully 
the quality of the gauge fixing and the importance of the Gribov copy problem, first on the 
propagator itself and then on the parameters that finally describe the functional form of the 
propagator. 

In a second direction, we investigate the potential influence that different definitions of the 
gauge potential in terms of the lattice link fields might have on the resulting propagator. 
One of the choices is strongly recommended by the explicit gauge invariance (transversality) 
of the emerging propagator. However, in the case when the longitudinal part does not vanish 
by construction, proper selection of the transverse component also leads to almost the same 
fit parameters. The coincidence becomes obviously better at larger f3. 

Going over from zero to finite temperature, it is important to realize that more struc- 
ture functions are necessary to the describe the finite-temperature case. One of them, the 
deconfinement-sensitive Dl, was studied already in Ref. Q, while the other, Dt, was found 
to be extremely sensitive to the Gribov problem. Only exerting extreme care can one expose 
the change going on at the deconfinement transition. 

The paper is organized as follows. In Sect. |H]we define the lattice model and the tensorial 
structure of the propagators, at T = and for T ^ 0. Sect. 11111 contains a discussion of the 
minimal Landau gauge which points out the topological aspects related to the gauge fixing, 
in particular the condition of having a minimal total length of Dirac strings. In Sect. II VI we 
report on the numerical algorithms for updating and gauge fixing which we have used in our 
investigation. The following Sect. |V|is devoted to an outline the results for T = 0. Sect. 
IVII describes the particular requirements of gauge fixing in the T^O case and contains our 
results for the propagators and Dt- Some conclusions are formulated in Sect. IVHI 

II. THE 3D COMPACT U(l) MODEL AND THE PHOTON PROPAGATOR 

For the compact U(l) model we chose the Wilson single-plaquette action: 



where 9 P is the U(l) field strength tensor represented by the plaquette curl of the compact 
link field Q\. The lattice is three-dimensional, and the basic degrees of freedom are the links 




(1) 



v 
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Ui = exp (i 9{). The measure of the link angles is flat over the interval —n < 9 < 71. The 
lattice coupling constant (3 is related to the lattice spacing a and the continuum coupling 
constant g% of the 3D theory, 

(3 = l/(agl). (2) 

Note that in three dimensional gauge theory the coupling constant g 3 has dimension mass 1 / 2 . 
Zero physical temperature is represented by symmetric lattices, L t = L s . 1 The lattice 
corresponding to finite temperature is asymmetric, L 2 s x L t , L t << L s . In the limit L s — > 00, 
the temporal extension of the lattice is related to the physical temperature, L t = l/(Ta). 
Using (0) the temperature is given in units of g\ in terms of the lattice parameters as follows: 



1 

g'i L t 



(3) 



Our simulations for zero temperature were performed mainly on a 32 3 lattice, those at finite 
temperature on a 32 2 x 8 lattice. 

The final discussion of the photon propagator will be given in lattice momentum space. 
Being always defined in a specified gauge, the propagator is written in terms of the Fourier 
transformed gauge potential, 

A S = - 1 Yexphm y K{ n » + 5<Wh A (4) 

which is a sum over a certain discrete set of points x = n + forming the support of 
Ag^ on the lattice. These are the midpoints of the links in /1 direction while n denotes the 
lattice sites (nodes) with integer Cartesian coordinates. The propagator is the gauge-fixed 
ensemble average of the following bilinear in A, 

D,M = {hAO ■ (5) 

Two identifications of were adopted in the literature and will be compared in our paper: 
the ang/e-definition 

A H+ i^ = On J (g 3 a) = log {U^) / (i £3 a) (6) 

and the sine-definition 

A n+ i^ = sin (6ft )IM ) I (g 3 a) = (U HtlM - JJtJ /(2 i g 3 a) . (7) 

The corresponding propagators will be denoted as D ^ 9 or D s ™, respectively. 

The lattice momenta p on the left hand side of (0) are related to the integer valued Fourier 
momenta k as follows: 

P M = -mn^, fc M = 0,±l,...,±^. (8) 
a l 



1 L t is the extension in the third (z) direction. 
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The lattice equivalent of p 2 = p 2 is in 3D 



^)^E(sinf) 2 . (9) 

At zero temperature, based on Euclidean rotational invariance, the continuum propagator 
would be expressible by functions of p 2 . The most general tensor structure is then the 
following one including two scalar functions of p 2 , 

d»M = p,M d(p 2 ) + (io) 

with the c?-dimensional (in our case d = 3) transverse projection operator 

Pr® = *i»- 2 f L - (11) 

The projector has the properties 

P,a(p) Pauif) = P^uip) , P,M=d-l. (12) 

The two structure functions D(p 2 ) and F(p 2 ) can be extracted by projection, on the lattice 
from Dp V {p) according to (J5J), as 

3 

and 

p 2 D ^ = d^T p ^ D ^ ■ (14) 

They are found, in the best case, only approximately rotationally invariant, i.e. individual 
momenta p might slightly differ in the function values D or F they provide, even if they 
have the same p 2 . Dense data points close together in p 2 might scatter rather than forming 
a smooth function of p 2 . 

In practice, using these definitions, we extract at first the function F(p 2 ), in the d = 3 
case through 

F(p 2 ) = plD n (p)+plD 22 (p)+plD 33 (p) 

+ 2p lP2 ReD 12 (p) + 2p lP3 ReD 13 (p) + 2p 2 p 3 ReD 23 (p) . (15) 

The imaginary parts of non-diagonal D^ v cancel in the sum and were omitted. Then, the 
function D{p 2 ) is obtained through 

F ) {p 2 ) = ^— — j- [(Ai + D 22 {p) + D 33 {p)) - F{p 2 )/p 2 ] . (16) 

If the Landau gauge would be exactly fulfilled, one would expect that F{p 2 ) = 0. On 
the lattice, in the case of the sine-definition for Ag^, this is actually the case as soon as 
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one of the Gribov copies is reached, with an accuracy which directly reflects the stopping 
precision of the gauge fixing procedure (as will be discussed below). In this case, a simplified 
definition in terms of the diagonal components D^{p) would be appropriate. In case of the 
ang/e-definition the structure function F(p 2 ) does not vanish, therefore all components of 
Dfj,u(p) contribute to D(p 2 ). 

For the finite temperature case, the propagator lacks 0(3) rotational symmetry. Now we 
have to consider two scalar functions, Dt and Dl instead of D, each multiplying the (d — 1)- 
dimensional transverse projection operator P T and the (d — l)-dimensional longitudinal 
projection operator P L , respectively. The scalar functions D T , D L and F depend now 
separately on the length of the space-like part of p with i = 1, . . . ,d — 1, 

P 2 =P? + ---+Pd-i, \P\ = V^, (17) 
and the "temporal" momentum p d : 

D,M = P*(p)D T (\p\,P d ) + P^P)D L (\ P \, Pd ) + • (18) 

The two projection operators are defined as follows (i,j — l,...,d—l): a transverse one, 



Pm = ^~ P -f> PJM = PJM = P? d (P) = o, (19) 

and a longitudinal one 

P.M = p»M - p»M ■ (20) 

Obviously, these projectors have the properties 

pUp) p^(p) = p»M > p Up) = d ~ 2 ( 21 ) 
PtM P «M = PlM , i%m = 1 (22) 

PlM p lM = °- (23) 

The scalar functions Dt and Dl can be extracted from D^ u (p) via 

PlM D»M = (d - 2) D T (\p\,p d ) (24) 

PlM d,m = D L (\p\, Pd ) 

= (d-l)D(\p\,p d )-(d-2)D T (\p\,p d ). (25) 

For d = 3 we can write down explicitly the definitions 

p 2 D T (\p\,p d ) = plD n (p) + p\D 22 {p) - 2 Pl p 2 ReD 12 (p) (26) 

and 

p 2 (p 2 +pt)D L (\p\,p d ) = p 2 {p 2 D 33 (p) - 2p lP3 ReD l3 (p) - 2p 2 p 3 ReD 23 (p)) 

+ V\ {p\DiM + p 2 2 D 22 {p) + 2p lP2 ReD 12 (p)) . (27) 

In the static limit, p 3 = 0, 

D L (\p\,p 3 = 0) = D 33 (\p\,p 3 = 0). (28) 
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This propagator for the case of the ang/e-definition and its fit in terms of mass and anomalous 
dimension as well as the changes at the deconfinement transition were discussed in Ref. 0. 

In the case of T = 0, the data to be presented below will be averaged over measurements 
of these quantities obtained for different k = (k%, &2, k^) giving rise to the same p 2 according 
to (JOJ). In the case T ^ we will show data averaged over different (fci,^) giving rise to 
the same p 2 according to ([17)1 . One should keep in mind that by this "trick" one enforces 
the rotational invariance "by hand", and the statistical errors are reduced as compared to 
measurements for individual ^-components of the propagator. 

In order to discuss the functional form of the propagator from the viewpoint of confine- 
ment effects, with confinement being induced by the monopole plasma, we will decompose 
the gauge fields into singular (monopole) and regular (photon) contributions on the level of 
the link angles, 

8*# = C + 9 ™ > ( 29 ) 

by a procedure to be described below. After the decomposition of the link angles is 
done, one may define the corresponding A p „ °\ _ and A™°Yl through the ang/e-definition 

(JSJ) or the sme-defmition (0), respectively. Then, by fast Fourier transform, the correspond- 
ing At hot and A™° no are evaluated. For each configuration and a certain set of momenta the 

k,fj, k,ii 

bilinears for the photon part ^ "t^pw ^ e m onopole part A™ ono A m ° no and the mixed bi- 

k,fi —k,u k,fi —k,u 

linear A phot A m ° no are formed. These are the observables which are associated - by averaging 

k,n ~k,u J 00 

over the Monte Carlo ensemble - to the propagators D p! l ot (p) (the photon or regular prop- 
agator), D™° no {p) (the singular propagator) and D™ xed {p) (the mixed propagator). These 
propagators are considered together with the full propagator which uses the original link 
angles 6^ before the splitting (j2Tlj) was performed. 

In the T = case, all these functions of p are then mapped by the projections (jl5H6j) 
to the scalar structure functions D phot (p 2 ), D mono (p 2 ), D mixed (p 2 ) and F phot {p 2 ), F mono (p 2 ), 
pmixed(p2y j n ^ e finite temperature case we proceed analogously. 

For the ang/e-definition of the vector potential we stress that the structure functions 
obviously satisfy exact additivity 

D(p) = D phot (p) + D mono (p) + 2 D mixed (p) (30) 
F(p) = F pho \p) + F mmo {p) + 2 F mixed {p) . (31) 



III. MONOPOLES, DIRAC STRINGS, AND THE MINIMAL LANDAU GAUGE 

In this section we partially follow Ref. [T3| • Gauge fixing to the Landau gauge means, for 
a given configuration Q\ = 9^^, to find a gauge transformations oj^ such that, for a gauge 
junctional 

T[Q) = Y,^i) (32) 
1 

the transformed gauge functional becomes maximal, 

mzxg(6,u), g(0,u) = F(9 {u,) ). (33) 
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Here vf ' denotes the gauge transformed gauge field 2 

9 -> 0H = 9 + duj 2 n = + diu + 2nk , k G Z , (34) 

where the integer number = fc(6j, u) for each link is chosen such that 9^ G (— tt, tt]. 

Instead of (J33J), following Ref. [lfj we use for the purpose of this section in the following 
the Villain form of the gauge condition 

minll^H 2 . (35) 

The Faddeev-Popov (FP) determinant is introduced by the following decomposition of 
unity 

1 = A FP [9;X\ [ Vue- Xlle( ^ 2 , (36) 



where A is the gauge fixing parameter. In order to achieve the gauge (J55)l we have to send 
A to infinity. In this case the limits of integration in (j3T?j) can be extended to ±oo since 
the saddle point approximation is exact in the limit A — > oo. Moreover, the integer valued 
variable k in (|34j) becomes effectively independent on 6 and us since the values of k for which 
q(v) ^ (—tt, tt] are exponentially suppressed. Therefore, in the limit of infinite A, the FP 
determinant (|36|) can be written as follows: 

oo 

A F l p [9;\] = f Vu expj-A||# + dcj + 27rA;|| 2 j. (37) 

oo fcex( ci ) 

Using the Hodge-de-Rahm transformation, k = 5A _1 dfc + dA _1 5/c, and making the shift 
oo — > uo — A~ l 5k we get 

oo 

A F ^[0;A] = const I Vu ^ exp|-A| \B + du + 2tt5A- 1 s\ | 2 | , (38) 

-OO »«( c 2) 
ds=0 

where we have changed the variables, s = dk. The integration over u gives 

App[6; A] = const exp|-A (dO + 2ir5A- 1 s, A^(d# + 2tt5A- 1 s) } . (39) 

ds=0 

To proceed further we separate the gauge field 9 into regular (photon) and singular 
(monopole) parts following Ref. 
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9 = 9 phot + 9 mono , 9 mono = 2ttA- 1 S P [j] , (40) 



2 Here and in the following we use the differential form notation on the lattice: (a, b) = ^ ; a;6/, ||#|| 2 = 
(9, 9). The operations d9 and 69 are the lattice curl and divergence, respectively. The Laplacian is denoted 
as A = Sd + dS. 
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where the dual one-form *j represents the monopoles on the dual lattice sites. The one- 
form on the dual lattice, p[j], defines the Dirac lines that connect the monopoles and anti- 
monopoles, S*p[j] = *j. 

The photon part phot is free of singularities while the monopole part Q mono contains the 
information about all monopole singularities: 

i-d[d^X = , ^d[d#™] 2w = J ■ (41) 



Here the DeGrand-Toussaint definition of the monopole |l2j was used. Substituting (|4f)j) 
into (f3HJ) we get, after a little algebra, 

A~^;A] = const exp{4A(j,A- 2 j)} £ exp{-\S gf (9 ph °\p[j] + s)} , (42) 



ds=0 



with 

S gf (6 phot ,p) = Uo phot + 2vrp, A^d^ * + 2tt{p[j] + s))\ . (43) 

The meaning of the last equations is the following. The gauge transformation (|3*4^) con- 
tains both regular (do;) and singular (k) parts. The former transforms the photon part of 
the gauge field while the latter changes the monopole part shifting the Dirac string (but 
leaving the monopoles j intact). We have already integrated out the regular gauge transfor- 
mations, therefore (|42j) depends explicitly on d6 phot which is invariant under regular gauge 
transformations, 6 phot -> 6 phot + duo. The sum in (J42)) over all possible shifts of the Dirac 
lines, 

*p\j]^*p\j]+*s, (44) 

corresponds to the integration over all singular gauge transformations (remember that *s 
is the closed line on the dual lattice, 5*s = 0). Thus (}4*2"|) is implicitly invariant under the 
singular gauge transformations as well. 

In the limit A — ► oo the only contribution to the FP determinant is given by the global 
minimum of the gauge fixing functional ()43|) with respect to the variations (|4^|) of the Dirac 
line, 

S™ n {6 ph °\j) = minS gf (e phot ,p[j] + s) . (45) 



ds=0 



If the photon field is absent, the minimum (|43j) is given by the Dirac line with minimal 
"Coulomb interaction" (c.f. ()43|) ). For a lattice monopole and anti-monopole separated 
along one axis this line is the shortest path connecting the pair. 

We substitute the FP unity ()36I42|) into the partition function of compact electrodynam- 



7T 

Z = Jv9e- s ^ ] (46) 
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then we transform the gauge field, 6> — > ^ and get the product of the gauge orbit volume, 
f Vuj, and the partition function within the fixed gauge, 

7T 

Zgf = J V6e- S ^- X ^ 2 A FP [8; A] • (47) 

— 7T 

Separating the gauge field on the monopole and photon parts as indicated by (jjQ)) and 
using the Hodge-de-Rahm transformation, one can show that 

| \9\ | 2 = (59 phot , Ar 1 56 phot ) + (&6 phot + 2vrp[j], A^d^' + 2vrj9[j])) - 4vr 2 (j, A~ 2 j) . (48) 

According to eqs. ()42I45I47|) the only non-vanishing contribution to the partition function 
in the limit A —>■ oo comes from the global minimum (}4*5|) in the gauge orbit. Comparison 
of (|45jl and (|4~H|) shows that this minimum is defined by the following conditions: 



together with 



s gf (e phot ,p[j}) = s^ n (e phot , P [j]) , (so) 



where S g f and S™p are given in (j4Hj) and (j43|) . respectively. 

In the continuum limit the condition (j4Uj) leads to the usual Landau gauge condition, 

d,Af ot = , (51) 

while the condition (joTIj) can be formulated as a requirement for the Dirac lines to form a 
configuration with as small as possible length: 

min length [p[j] ) . (52) 
S*p\j]=Sj V / 

Indeed, the Dirac lines, *p[j] = *Pi[j] + *P2[j] + • • • , correspond to the singular 5-functions 
in the continuum limit. Here *pi correspond to mutually unconnected pieces of these lines. 
The self-interaction of the Dirac lines in ()50|) . Y2ii*Pi[j]i ^ 1 *Pi[j])i contains the term 
a length (with a logarithmically divergent coefficient a) plus finite terms. The 

"Coulomb interaction" of different pieces of the Dirac lines, [*Pi[j], A~ 1 *pk[j]), i ^ k, as well 
as the contribution of the regular fields into the condition (jBHj) are finite in the continuum 
limit. Thus the only essential contribution to the condition (jBUjl in the continuum limit is 

given by the term a J2i length = a length (p[j]J , which gives the condition (p)2"j) . 

Thus we conclude that in the continuum limit the minimal Landau gauge for the compact 
gauge fields is reduced to the local gauge condition (foT|) for the regular fields and a non-local 
condition (|52|) - the requirement for the total length of the Dirac lines to be as small as 
possible - for the singular fields. This result can easily be generalized to the 4D case. 
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IV. NUMERICAL ALGORITHMS USED IN THE ANALYSIS 
A. Monte Carlo Updating 

The Monte Carlo algorithm in use for this investigation is a mixture of local and global 
updates. The local Monte Carlo algorithm is based on a 5 hit Metropolis update sweep in 
an even-odd fashion, alternating with a microcanonical sweep, also in checkerboard mode. 
Both together are considered as one local update. After three local updates the Metropolis 
step width is eventually tuned to keep an acceptance in the range between 40 % and 60 %. 

For better ergodicity, in particular in the presence of an external field (considered in 
Ref. 0), global updates have also been included, following the ideas of Ref. [l3| . In the 
equilibrium regime, after every three complete local updates, a global refreshment step is 
attempted. We try to add one unit of flux to the dynamical gauge field, with random sign 
in one of the three directions randomly selected. The proposed flux addition is subject to a 
global Metropolis acceptance check. 

For example, one unit of flux in the pv plane is introduced with the help of the following 
gauge field shift [13] 6^ -> [9 S}fl + 6^} mod 2vr : 

7T 

9x,v = y-(2 — — 1) , 3jV = for x u ^ L u , 
2-n 

9s,n = —L u (l-x u ), 9 3>P = Q, p^fi,u. (53) 
m 

The acceptance rate of the global step changes with /3 in a different way, depending on 
the lattice geometry. In our T = studies (on 32 3 lattices) we found that the acceptance of 
global update steps drops (more rapidly than exponentially) from 0.48 at (3 = 1.0 to 0.0056 
at (3 — 2.0. For higher (3 essentially no global offers are successful. 

In our T ^ studies on 32 2 x 8 lattices, however, we found the acceptance changing 
smoothly (nearly exponentially, even across the deconfining transition) from 0.58 at (3 = 1.0 
to 0.18 at (3 = 3.0. A closer look reveals that the higher acceptance rate is due to more 
frequent global changes of the flux penetrating the 12 plane (i.e. magnetic flux direction). 

In summary, one total Monte Carlo update cycle consists of three cycles of local update, 
each consisting of a Metropolis sweep followed by a microcanonical sweep, interchanging 
with a global update as described above. 

In the finite temperature case, the measurement of Dt turns out to be highly sensitive 
with respect to the insufficient removal (by the gauge fixing procedure and its repetitions, 
see the next subsection) of Dirac strings wrapping around the third direction. This is a case 
where the results with and without global updates, mainly adding and subtracting fluxes 
through the 12 plane, differ. We comment on this problem and how to deal with it in 
Section ED 



B. Landau Gauge Fixing 

The Landau gauge was chosen first of all because it is the most popular gauge to define 
a gauge field propagator. In this gauge the gauge propagator (jSJ) is expected to satisfy the 
transversality condition 

F(q 2 ) = q, D, v {q) q v = . (54) 
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In the case of zero temperature, for example, this allows to describe the propagator by a 
single function D(q 2 ) alone, defined by 

D»M = (V - ^) D(q 2 ) . (55) 

We will see that for any practical implementation of the Landau gauge (p>fj) is slightly 
violated. This degree of violation can, however, be easily controlled by sharpening the 
convergence criteria of the gauge-fixing algorithm. More important is the remark that in 
the case of the ang/e-definition of the vector potential (JHJ) it is really necessary to select the 
transverse part by projection using (fTB} and (fTB|) . Then it is interesting to see where (e.g. in 
momentum space) the violation of transversality, quantified by the longitudinal propagator 
F, is coming from. 

There is a second reason to choose the Landau gauge. We intend to split the gauge field 
into a regular (photon) and a singular (monopole) part by reconstructing the field due to the 
Dirac plaquettes which, on the other hand, are forming the monopoles. This reconstruction 
becomes unique in the Landau gauge. 

In order to implement the gauge fixing condition (IHHj) we have chosen a mixture of over- 
relaxation and non-periodic gauge transformations |l4j . both applied in alternating order. 

Iterative overrelaxation has to be practiced in a checkerboard fashion. Starting, say, with 
the odd sublattice we have first to find for each odd site n a suitable which maximizes 
the following function of uj 

G l S c (6, u) = ( cos (^ - + cos(0a-A** + un)) , (56) 

which represents the part of Q actually depending on cj^, and second to perform immediately 
the updatings of the neighboring link angles 

This can be done simultaneously for half of the sites, namely n 6 A ot ^. Afterward the same 
procedure is applied to the even sublattice. One odd/even pair of gauge updates constitutes 
one single iterative overrelaxation step. Each overrelaxation iteration is followed by a zero 
mode subtraction (to be explained in the next subsection). 
The angle lo^ can be easily found as 

£jsin(0^) -sin(^_ A/ J 

tan(u%) = j ; 

Z)J cos(0jj,/J + cos(6> 

These gauge angles are multiplied by the overrelaxation factor, — > r/uj^ and bounded by 
\ujfi\ < it before the iteration (|H7j) is performed on all links. A good overrelaxation parameter 
was found to be r\ = 1.8. According to our experience from studies on 16 3 lattices this r] 
leads to fastest convergence, almost independently of /3. This value has then been applied 
for all iterative gauge fixings. 

The overrelaxation will usually be stopped if in the last overrelaxation step the average 
increase of the gauge functional T [Eq.(j32J)] per link is found to be less than 10~ 6 . After 



(57) 



(58) 
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this was discovered, the gauge fixing procedure ends with a final zero mode subtraction (see 
below). 

At any local extremum of T [Eq. (j32jl ] the following condition would be satisfied everywhere 
on the even and odd sublattice: 

(W)n = £(4t4-iA, - A n _^) = -L ^(sin(^) - sin(^)) = . (59) 

|U A 4 

Having the vector potential localized on the midpoints of links, its divergence is naturally 
defined on sites n. Exact vanishing of the divergence of can be expected in the result of 
Landau gauge fixing only for the sine-definition of the vector potential Ap. 

The algorithm outlined above will not in general lead to the absolute (global) maximum 
of the gauge functional T ()32|). Typically it will get stuck in one of the local maxima, which 
are called Gribov copies of the true maximum. This is the so-called Gribov problem. It is 
partially cured by repeating the same gauge fixing procedure, applying it to random gauge 
copies of the original Monte Carlo configuration 9 MC , assuming that one of these copies 
might be situated in the basin of attraction of the true maximum. The number of gauge 
equivalent configurations produced to restart the gauge fixing is denoted as N G , and the 
iterative gauge fixing generically leads to really different maxima. We have then to content 
with the best out of all N G + 1 local maxima of the gauge functional. The convergence, with 
increasing Nq, of a gauge dependent quantity evaluated on a given gauge field ensemble with 
the help of the best Gribov copy gives an indication of the sensitivity of this quantity with 
respect to the misidentification of the true maximum. We have applied this philosophy to 
two sets of data, the propagator data themselves at large or small momenta, and to the fit 
parameter emerging from a fit of the gauge boson propagator. One should not be surprised 
that the number Nq, which is necessary to achieve uniform convergence of the propagator 
in momentum space and/or of the fit parameters, differs strongly between zero and finite 
temperature. At finite temperature there are strong differences between and D?- 

In section 11111 we have stressed the importance to find a local maximum of (|32|) accom- 
panied by a minimal length of Dirac strings. Within our implementation, the start from a 
new random gauge copy is done in the hope in producing a new Gribov copy reachable from 
the previous ones only by a discrete gauge transformation. We have monitored the number 
of Dirac strings in each of the local maxima of the gauge functional T (|32p. Each time the 
recent best value of the gauge functional was replaced by a better (higher) one, the number 
N D of Dirac strings detected in the correspondingly best Gribov copy decreased compared 
to the previous one. 

In Fig. Q we show scatter plots relating the gauge functional T to the number of Dirac 
strings per direction Nf (i = 1, .., 3)) for two thinned ensembles of 100 out of 1000 equilib- 
rium configurations on 32 2 x 8 lattices, each of them entering with 101 Gribov copies (local 
maxima of J 7 ). Fig. E^a) refers to /3 — 2.0 (confined phase) and shows all local maxima 
(Gribov copies). The isotropy of Dirac strings is clearly visible. Fig. ^b) shows how in the 
deconfined phase the number of temporal Dirac strings is correlated with the locally 
maximal value of T . Fig.^c) shows it for the Nf and (number of spatial Dirac strings). 
In both cases a slight tendency to clustering near multiples of 8 or 32 (the respective peri- 
odicity of the lattice) is visible, in particular for the highest values of T . The selection of 
the best copies restricts the ensemble to bigger T and smaller Nf as shown in Fig. C^d). In 
the confinement phase the Dirac strings are still approximately isotropic. In the deconfined 
phase, however, the highest values of T are correlated with low multiplicity of N± and N®, 
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FIG. 1: Scatter plots of the average gauge functional J-{6) per link and the related number of Dirac 
strings in the temporal (iVjP) and/or the two spatial (N^ and N^) directions for 100 configurations 
on a lattice 32 2 x 8: (a) Nj^ , Nf and vs T(9) in the confinement phase (f3 = 2.0), each 
configuration represented by 101 Gribov copies; (b) vs J-(6) in the deconfinement phase (/3 = 
2.6); (c) same as in (b) for and ; (d) the two samples represented only by the best out of 
101 Gribov copies. 



and corresponding to zero, with one or a few periodically winding Dirac strings more 
frequently in the temporal direction. 

We stress, that the configurations of which the gauge fixing was investigated here in some 
detail, had been produced with the update algorithm including global updates. Winding 
Dirac strings of some life time are also produced in ensembles that are generated without 
global updates, but less frequently. Therefore, the inefficiency of the random gauge trans- 
formation in exploring more of the gauge orbit is a handicap also if global updates are 
suppressed. We will see later that certain problems which show up in the finite temperature 
propagator D T can be ameliorated, but not completely cured, by abandoning global updates 
and increasing the number of Gribov copies. 

We had not the opportunity to reconstruct exactly the spatial conformation of the Dirac 
strings in each gauge copy. But for each gauge copy we can determine whether multiple 
Dirac strings (running along a certain direction with same or different orientation) can be 
definitely excluded or not. At low (3 this can never be excluded, but at high (3 this could be 
used as an additional criterion to reject Gribov copies which contain double Dirac strings 
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(in case of two oppositely oriented strings). 

Although the stopping criterion was formulated for the global increase of T we did not 
find a systematic local variation, depending on the distance from a monopole, of the violation 
of (|59|). expressed by the quantity 



v v 



— A 



(60) 



This suggests that the differential gauge condition (|59)l is uniformly approaching zero. 

In contrast to this, we found a systematic local variation of the local gauge functional 
itself, 



(61) 



with the distance from a monopole. Near to a monopole, its value is suppressed compared 
to the bulk average. This is illustrated in Fig. |2]for (3 = 1.0 and (3 = 2.0. 




FIG. 2: Suppression of the local gauge functional T l £ c near to a monopole for j3 = 1.0 and j3 = 2.0. 
For better presentation part of the data points are not plotted. Notice the f3 dependence of the 
bulk average of 



C. Zero Mode Subtraction 

There are certain modes of the gauge field which are not suppressed by the action. If 
they are not taken care of properly, the gauge field propagator is known to be spoiled [14j |. 
For instance, adding some constant to all link angles 0^ i(U in one particular direction does not 
change the action (0), as neither does adding a multiple of 2tt to one particular link angle. 
Non-periodic gauge transformations can be considered to implement these changes. We have 
used it in our simulations to immediately eliminate zero modes related to the appearance of 
the volume average of the link angles in a given direction 



n 



(62) 



15 



where |A| is the lattice volume. We subtract this volume average from each link angle 9^ 
in fi direction, 

9n,^9n^-9; + 27rk, (63) 

with integer k chosen such that 9 f i G (— 7T, tt]. This operation can be imagined as resulting 

from a gauge function gn — — 9^n^ (no summation!) which is in general non-periodic. This 
special sort of gauge fixing is completed when zero modes in all d directions were subtracted. 

The zero-mode subtraction step is applied after each overrelaxation step. Thus it is 
always performed before the measurements are done on the gauge-fixed configuration. 

In our runs we measured the photon propagator after each 10th total Monte Carlo update 
cycle to avoid autocorrelations as much as possible. Typically we used 500 gauge-fixed 
configurations per data point for the 32 3 lattice and about 2000 configurations for the 32 2 x 8 
lattice. 



V. THE ZERO-TEMPERATURE PROPAGATOR IN LANDAU GAUGE 

In coordinate space, the gauge boson propagator that is studied, is read 

(™ + \fl - ™ - \vj = (A A+ i^A H+ i PjU ) (64) 

with the angle- or sine-definition of A. For brevity, we will refer to this propagator later 
in momentum space as F) ^ 9 or D 8 ™. The two scalar functions occurring in momentum 
space (|10|) we will denote as D ang and F an9 or D sm and F sm , respectively. The last one, 
F sm should vanish in the Landau gauge. We have observed that this is indeed the case 
with an accuracy determined by the stopping criterion of the iterative overrelaxation. For 
the sine-propagator the transverse part could be calculated directly just by evaluating and 
summing the diagonal components appearing in (|16p. When the longitudinal propagator 
vanishes only approximately, one can extract the transverse propagator following ()151 [TBI) . 

In Fig.|3]we show the different forms of the transverse propagator D sin and its components 
as well as the vanishing longitudinal propagator F sm . These data were obtained on a 32 3 
lattice for (a) j3 = 1.0 and (b) (3 — 1.8, with N G = 20 Gribov copies evaluated in addition 
to each original Monte Carlo configuration. The data at other /3 were produced under the 
same conditions. In Fig.0]the same is presented for D an9 , its components and non-vanishing 

pang 

For the ang/e-definition of A^, the decomposition into components is strictly additive 
(|31|). We have observed that the longitudinal propagator F ang and its components are well 
described by a form F(p 2 ) = P p 2 , where P is a constant. We find that F an 9^ hot essentially 
coincides with the full F an9 . The size of F an9 and its photon component is of the same 
order of magnitude as the transverse part at (3 = 1.0, while the monopole and mixed part 
are one order of magnitude smaller. At (3 = 1.8, the size of F ang and its photon component 
is an order of magnitude smaller than at (3 = 1.0, while the monopole and mixed part are 
negligible. 

For the sine-definition of A^, the monopole part of the transverse propagator D sm < mono 
has a maximum in the low momentum region (which moves more and more towards p 2 = 
with higher (3) before it drops towards p 2 = 0. For the ang/e-defmition, a maximum of 
D ang,mono develops only for (3 > 1.0. 
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FIG. 3: The Landau gauge sine-propagator as a function of p 2 measured on a 32 3 lattice : for 
the transverse part, p 2 D sm ((a) at f3 = 1.0 and (b) at f3 = 1.8) we show the full propagator, the 
singular (mono), the regular (phot) and the mixed contribution for comparison. In addition, we 
show the (vanishing) longitudinal propagator F sm . The data represent the evaluation of Nq = 20 
gauge copies. For better presentation part of the data points are not plotted. 
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FIG. 4: The same as in Fig. |3] for the Landau gauge angle propagator. Notice the nonvanishing 
longitudinal propagator F ang decreasing with growing (3. 



Summarizing, we have observed that the longitudinal part F is either zero (for the sin 9 
definition of the gauge field in the correlator) or it is non-zero, and then it coincides with 
its photon part (for the 9 definition). Therefore, at zero temperature the Landau gauge 
propagator D a ^f is not completely transverse. This is entirely due to the difference between 
the definitions of the vector potential. This discrepancy, expressed by the nonvanishing 
pang ^ k ecomes ameliorated at higher (3. 

For both definitions of we find that the regular (photon) part of the transverse prop- 
agator is singular at p 2 — > like D phot ~ 1/p 2 , while the full transverse propagator is not. 
Following Ref. 0] we try to describe the two by functions of the form 

Z m 2a 
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and 



D phot {p 2 



Z p, 



hot 



(3 p 



+ C ph 



ot 



(66) 



The model function (f63J) is similar to some of the functions discussed in Refs. [lol flfij ] where 
the propagator in gluodynamics was studied. In the case of T = we expect these two 
curves to differ at all (3 in order to accommodate the (permanent) confinement property of 
the model. 
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FIG. 5: The best fit parameters for the zero temperature sine- and cmg/e-propagators as functions 
of (3: (a) the anomalous dimension a, (b) the mass parameter m with the theoretical prediction 
(|BTj) as dashed line, and (c) the parameters Z and Z p h ot (using ()65I66|0 : (d) the total monopole 
density and the monopole density in charged clusters with the theoretical prediction 



Fig. El^a) shows the anomalous dimensions a first increasing in the low-/? region. The 
anomalous dimension for the angle- and for the sine-propagator behave quite similar to 
each other (the dimension for the ang/e-propagator is a bit smaller). As (3 gets larger 
than (3 ~ 1.5 the anomalous dimensions start to descent towards zero. This indicates that 
the anomalous dimension is not only a function of the monopole density which decreases 
monotonically with growing (3 for all values of the coupling. The (cluster) structure of the 
monopole configurations may play a significant role for a. 

In Fig.^b) the mass parameters m are presented for the angle- and the sine-propagator 
according to as function of (3. Both masses are almost equivalent. For the mass there 
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exists a theoretical prediction due to Polyakov [1 




(0)}, 



(67) 



where (3y is the Villain coupling constant 



Pv(/3)= 2 log 



( 



h((3) 



) 



(68) 



Jo(/3) and Ii(/3) are the modified Bessel functions and (3 is the Wilson action coupling 
constant appearing in (JTJ). The prediction (|H7jl is valid for a dilute monopole gas. The 
agreement between the two data sets and the theoretical curve is very good. The small 
deviation at lowest (3 can be attributed to the violation of the dilute gas approximation. 

Fig. Efc) shows Z(/3) and Z phot (j3) for the two definitions. We observe that Z tends to 
unity at large /3, whereas Z phot {(3) « 1 for all (3. The strong deviation of Z from unity at 
small j3 can be interpreted as a field renormalization by monopoles. 

The simplest quantity characterizing the monopoles is the monopole density, p = ^Tj C j c / 
|A|, where the monopole charge j is defined in Eq. (14 lj) . Note, however, that a general 
monopole ensemble may contain lattice artifacts which at zero temperature are realized in 
the form of ultraviolet monopole-antimonopole pairs. Following Ref. |5| we remove these 
lattice artifacts using a cluster analysis. For our purposes, clusters are defined as connected 
groups of monopoles and anti-monopoles, where each object is separated from at least 
one neighbor belonging to the same cluster by a distance less or equal to -R max - We use 
-^max = 3 a 2 which means that neighboring monopole cubes should share at least one single 
corner. The cluster is called charged if the total charge of its constituent monopoles is 
non-zero. This includes isolated monopoles and anti-monopoles. 

In Fig. E(d) we plot the measured total monopole density and the density of monopoles 
residing in charged (physical) clusters. The charged fraction is well described by the theo- 
retical formula for the monopole density, 



which is a lattice version |5[ of the Polyakov formula [l| . According to Fig. \5\ both the 
monopole density, Debye mass and the deviation of the coupling Z from unity are descending 
functions vanishing in the limit /3 — > oo. 

The contact terms contained in D sin and D ang are not shown here. The photon part of 
both D sm and D ang vanishes perfectly. The full propagator in both cases contains contact 
terms, C sm and C an9 , which deviate from zero for smaller (3, whereas always C sm ((3) <C 
C an9 {(3). 

The discussion of Figs. El and |U and of the (3 dependence of the fit parameters in Fig. E] 
was based on the zero-temperature propagators (and its components) obtained throughout 
with Nq = 20 Gribov gauge copies. The dependence on the number of gauge copies was 
investigated carefully for the case of the angZe-propagator. In the result of this study, the 
default choice of N G = 20 for the gauge-fixing procedure at T = was established. 

In Fig. El we show different aspects of the approach to the N G — > oo limit, for (3 = 2.0 
as an example. Fig. a) demonstrates that the transverse propagator evaluated with 10 or 
20 gauge copies is essentially the same, but that the naive evaluation (with N G = 0) would 
clearly overestimate the propagator over the whole momentum range. Fig. EJb) shows this 




(69) 
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FIG. 6: The dependence of the zero-temperature ang/e-propagator on Nq for a 32 3 lattice (at 
f3 = 2.0 as an example) : (a) the full transverse propagator for Nq = 0, 10, 20 in the full momentum 
region; (b) its behaviour as function of Nq at five selected momenta; (c) the same for the singular 
(mono) part of the transverse propagator; (d) the same for the regular (phot) contribution to the 
transverse propagator. Again for better presentation part of the data points are not plotted. 



in more detail for five selected momenta. It becomes clear that the dependence is strongest 
in the region of small momenta, in particular the region below the peak. The dependence is 
strong for the singular part presented in Fig.H^c). Again, there is almost no change between 
Nq = 10 and 20. As can be seen from Fig. Efd), there is almost no Nq dependence in the 
photon part of the transverse propagator. 

We present the resulting dependence of the fit parameters on Nq in Fig. [7[ again for 
(3 = 2.0. The anomalous dimension a shown in Fig. Ufa) drops within 10 % which indicates 
that the anomalous dimension is sensitive with respect to the minimization of the Dirac 
strings that is achieved by better and better Gribov copies. The mass presented in Fig. 0T>) 
does not change with Nq which confirms that it is mainly determined by the monopole mass 
(i.e. density). 3 Fig.Efc) shows Z for the full transverse propagator and Z phot for the photon 
part. It is not surprising that the parameter associated with the photon part do not change. 
The parameter describing the full transverse propagator decreases, again by m 10% within 
Nq < 20. One can also observe a slight dependence of Z on Nq for Nq > 20. However, the 



3 Let us recall that the monopole positions in a given configuration are gauge independent. 
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FIG. 7: The best fit parameters for D an9 at zero temperature as a function of Nq (for /? = 2.0 as 
an example): (a) anomalous dimension a, (b) mass parameter m and (c) renormalization constants 
Z for D an 9 and D an9 ' phot . 



dependence is indeed very small (about 1%) and not essential for our qualitative discussion. 

From the comparison of the sm^propagator and the ang/e-propagator at zero tempera- 
ture we conclude that the fits of the transverse part give more or less the same parameters, 
with a j3 dependence (if appropriate) which is in accordance with the monopole density. 
Strict transversality itself is guaranteed only in the case of the sine-propagator. In the 
case of the ang/e-definition an appropriate transverse part has to be extracted by projection 

flnnEp. 



VI. THE FINITE-TEMPERATURE PROPAGATOR IN LANDAU GAUGE 

In this Section we report our investigation of the properties of the gauge boson propagator 
at finite temperatures. With respect to the distinguished direction fi = 3, the propagator 
D can be separated into transverse and longitudinal components (see Section |n] for the 
definitions) denoted as and Dl, respectively. We are working on the lattice 32 2 x 8, in 
line with Ref. Q, where only the Dl component of the ang/e-type propagator (with p 3 = 0, 
Nq = 20 and a limited statistics of 500 measurements) was studied. 

The transverse component of the propagator, Dt, describes the spatial degrees of free- 
dom while the longitudinal component, Dl, contains gauge fields in Landau gauge in both 
temporal and spatial directions. The finite temperature propagator data are analyzed again 
for p 3 = 0, as a function of p 2 . In that case Dl is constructed only from temporal degrees 
of freedom which, in particular, are responsible for the confinement phenomena. We have 
fitted the data for both components of the propagator using the fit function (p3o"|) invented 
first in Ref. to describe D 33 . 

First we repeated the investigation of the Gribov copy dependence of the propagator 
components Dl and Dt, this time for the sme-definition of the propagator, similar to that 
conducted for the zero-temperature case with the ang/e-definition. The results are summa- 
rized in Fig. |S1 at (3 values near (below and above) the phase transition. After a few Gribov 
copy attempts the longitudinal component is almost insensitive to the number Nq. This 
can be seen from the left panel of Fig. |H] The fitting parameters o>l, m^, Z L and Cl (the 
latter is not shown here) are rapidly converging and become almost independent of Nq for 
Nq > 7. The results at large number of Gribov copies are not sensitive to the fact whether 
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FIG. 8: The Gribov copy dependence of the best fit parameters for the Dl and Dt propagators: 
the anomalous dimension a (a,b), the mass parameter m (c,d) and the parameters Z (e,f). The 
label "local" corresponds to measurements in Monte Carlo cycles without global updates and with 
Nq = 100 extra Gribov copies. 



we have suppressed global updates (we made only local updates, see the label "local ") or 
not. 

The transverse component, however, is strongly dependent on the number of Gribov 
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copies as it can be seen from Fig. |H| All fit parameters «t, my and Zt are descending 
functions of the number of Gribov copies N G with global updates included. At N G = 100 
the plateau is not yet reached. Moreover, the results are sensitive to whether or not global 
updates are included, in particular at high values of [3. The measurements with only local 
updates lead to significantly lower fit results. In deconfinement we would expect vanishing 
rriT and aj), and Zt — > 1. 

The reason for that behavior might be explained as follows: On one side, the "best" gauge 
functional is realized in gauge-fixed configurations without Dirac lines wrapping around the 
lattice (predominantly in the short temporal direction). Such Dirac lines are continuously 
created and destroyed by the Monte Carlo process, even if global updates are not attempted. 
The level of "noise" due to wrapping Dirac lines is higher if global updates are included in 
the Monte Carlo process, which in general would improve the ergodicity of the system, 
but it represents a problem also if only local updates are used. The presence of wrapping 
Dirac strings mimics a finite /^-independent lattice mass tjit at larger (3, the value of which 
decreases only with increasing temporal extent Lt- So in the limit of vanishing lattice 
spacing the dimensionful mass would diverge. 

This "Dirac noise" represents a serious challenge for the gauge fixing algorithm. The 
deterministic part (overrelaxing steepest descent method) described in Section TlV Bl cannot 
remove it. Unbiased random gauge transformations applied to get new start configurations 
for the deterministic search for further Gribov copies are obviously not effective enough 
to reduce the "Dirac noise". A simulated annealing Monte Carlo series of random gauge 
transformations with the total length of Dirac strings as "gauge action" [l7| seems to be 
more appropriate for selecting new start configurations for the final steepest descent search. 

Having these difficulties in mind, we decided to use in the final measurements at finite 
temperature, for j3 = 2.0 and larger, only local updates before gauge fixing and to perform 
Ng = 100 Gribov copy attempts. For both sine- and ang/e-propagator measurements deep 
in the confinement phase (below (3 = 2.0) we used N G = 20 and global updates where the 
results of the zero temperature analysis for the Gribov copy dependence is applicable and 
the fit parameters of both Dl and Dt have to agree within accuracy and should be similar 
to those for the T = transverse propagator D. Nevertheless, we have to admit that the 
results for the transverse propagator D T should be understood only qualitatively. 

The results for best fit parameters for and Dt are presented in Fig. El Let us begin 
with D^. At the critical point both the anomalous dimension and the mass vanish 
while the renormalization parameter Zt meets the corresponding parameter for the pertur- 
bative photon, Z^ hot . This behaviour is characteristic for both angle and sine types of the 
propagator and extends our results in Ref. jjfl. Note that also here the mass parameters for 
sine- and cmg/e-propagators coincide with each other. However, the anomalous dimensions 
for these cases differ slightly from each other while the renormalization parameters are sig- 
nificantly different, and the renormalization factor for the ang/e-propagator is bigger than 
that for the sme-propagator. The latter is expected because | sin#| ^ \9\. 

The corresponding quantities for Dt behave differently for the angle- and sme-definitions 
of A^, with the remarkable exception of the mass parameter. For example, the anomalous 
dimension for the ang/e-propagator vanishes in the vicinity of the critical point and 
beyond, while the same quantity for the sine-propagator does not vanish. We explain this 
behaviour as due to insufficient gauge fixing as it can also be guessed from our previous 
analysis. The same reason explains the fact that the masses uit for both definitions of the 
photon propagators - being remarkably similar - are not vanishing at the critical point. 
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FIG. 9: The best fit parameters for the non-zero temperature sine- and ang/e-propagators as 
functions of (3: the anomalous dimension a (a,b), the mass parameter m (c,d), the parameters Z 
and Z p hot (e,f)- The left column corresponds to D^, the right to Dt- The fits of full propagators 
are done with the help of the photon contribution is fitted by (|Bd]) . 



Finally, for both definitions of the propagators the renormalization constants do not 
approach the corresponding Z T hot at the critical point. In order to get a reliable behaviour 
of Dt part of the propagator one should drastically increase the number of Gribov copies 
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used in the gauge fixing. For the time being this is beyond our computing capabilities. The 
situation could be improved using a variant of the mentioned simulated annealing Monte 
Carlo series of random gauge transformations in order to choose more appropriate initial 
gauge transformed configurations before fixing the gauge. 

VII. CONCLUSIONS 

We have studied the gauge boson propagator in cQED 3 both at zero and non-zero tem- 
peratures. We have found that the propagators in all cases under investigation can be fitted 
by which is the sum of the massive propagator with an anomalous dimension plus a 
contact term. Similarly to the case of Dl at finite-temperature 0], a nonvanishing anoma- 
lous dimension a is also found at T = 0. Moreover, the fact of existence of the anomalous 
dimension is not associated with a particular type of gauge boson propagator. We have 
studied angle- and sine-types of propagators and the corresponding anomalous dimensions 
are nonvanishing and have a similar behaviour as the functions of f3. 

The existence of the anomalous dimension depends on the presence of the monopole 
plasma, but it is not directly proportional to the monopole density below /? = 1.5. In the 
confinement phase the monopole plasma is present at any coupling of the system and the 
density of monopoles is a monotonously decreasing function of the lattice coupling (3. A 
similar behaviour is observed for the anomalous dimension in the case of the Dl and D T 
propagators. The dimension extracted from the Dl component of the propagator is 
vanishing in the vicinity of the phase transition for both definitions {angle and sine) of the 
propagator. However, this does not happen for the sine-definition of the Dt propagator. 
We associate this result with insufficient number of the Gribov copies used in the gauge 
fixing. The Dt propagator requires many more Gribov copies for the gauge fixing than does 
the D L propagator. 

Concerning the other parameters of the fits, the mass extracted from the propagator at 
zero temperature and from the Dl propagator at non-zero temperature does not depend 
on the definition of the propagator. The mass for the T = case is perfectly described by 
the Polyakov formula. The tul mass vanishes at the phase transition point as was expected 
from the disappearance of the monopole plasma at the critical temperature. Beyond the 
phase transition point, the mx mass measured in this paper does not behave in a physical 
way due to the severe Gribov copy problem. 

Finally let us comment on the continuum limit of the measured quantities. The continuum 
limit of the cQED 3 corresponds to (3 = 1/(<?J a) — ► oo holding the dimensionful gauge 
coupling g 3 fixed. According to Polyakov, non-perturbative quantities such as the Debye 
mass and string tension can be expressed in terms of g% and the monopole density p, which 
generally might be independent quantities. However, this is not true for compact U(l) where 
both g 3 and p depend on a single parameter, the lattice coupling (5. Therefore, in the limit of 
vanishing lattice spacing the monopole density and other non-perturbative quantities such 
as the Debye mass, string tension and anomalous dimension also vanish exponentially as 
~ exp{— const (3} (cf. Eqs. ()6 7169(0 . However, in more realistic models (like the Georgi- 
Glashow model) the monopole density and the lattice spacings are indeed independent and 
the monopole density should survive in the continuum limit. According to our results, this 
implies that in the continuum limit of such theories a non-zero anomalous dimension in the 
photon propagator can be expected. 
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